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Abstract 

In this paper, we introduce a Markov process whose unique invariant 
distribution is the Curie-Weiss model of self-organized criticality (SOC) 
we designed and studied in [4]. In the Gaussian case, we prove rigorously 
that it is a dynamical model of SOC: the fluctuations of the sum Sn( •) 
of the process evolve in a time scale of order i/n and in a space scale of 
order and the limiting process is the solution of a "critical" stochastic 
differential equation. 
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1 Introduction 


In [4] and [10], we introduced a Curie-Weiss model of self-organized criticality 
(SOC): we transformed the distribution associated to the generalized Ising Curie- 
Weiss model by implementing an automatic control of the inverse temperature 
which forces the model to evolve towards a critical state. It is the model given 
by an infinite triangular array of real-valued random variables (^^)i<fe<n such 
that, for all n > 1, {X^,..., X^) has the distribution 

^ xl + --- + xl ) 11 Wi), 

where p is a probability measure on R which is not the Dirac mass at 0, and 
where is the normalization constant. We extended the study of this model 
in [11], [12] and [13]. For symmetric distributions satisfying some exponential 
moment condition, we proved that the sum Sn of the random variables behaves 
as in the typical critical generalized Ising Curie-Weiss model: the fluctuations 
are of order and the limiting law is C exp(—Ax^) dx where C and A are 
suitable positive constants. Moreover, by construction, the model does not 
depend on any external parameter. That is why we can conclude it exhibits the 
phenomenon of self-organized criticality (SOC). Our motivations for studying 
such a model are detailed in [4]. 

This model describes interacting elements in thermodynamic equilibrium. How¬ 
ever self-organized criticality seems to be a dynamical phenomenon, as is high¬ 
lighted by the archetype of SOC : the sandpile model introduced by Per Bak, 
Chao Tang and Kurt Wiesenfeld in their seminal 1987 paper [1]. That is why, 
in this paper, we try to design a dynamical Curie-Weiss model of SOC. 

We choose to build a dynamical model as a Markov process whose unique invari¬ 
ant distribution is the law of (a modified version of) the Curie-Weiss model of 
SOC. One way of building such a process is to consider the associated Langevin 
diffusion (see [16] for example). 

The model. Let p be a function from R to R which is even and such that 
the function exp(2(p) is integrable over R. We suppose that there exists C > 0 
such that 

Vx G R xip'{x) < C{\ + x^). 

We denote by p the probability measure with density 

XI—>■ exp(2(p(x)) exp{2(p{t)) dt^ 


with respect to the Lebesgue measure on R. We consider an infinite triangular 
array of stochastic processes (X^(t), t > 0)i<fc<„ such that, for all n > 1, 


is the unique solution of the system of stochastic differential equations: 


dXi{t) = p'{Xi{t))dt + dB,{t) 

1 ( Sr^jt) 
2\Tn{t) + l 


xm 


Suit) 
Tn{t) + 1 


^ dt 

j G {l,...,n}. 


(S)?) 
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where (i?i,. •., Bn) is a standard n-dimensional Brownian motion and 

vt>0 Snit)=xlit)+---+x::it), Tnit) = {x^nit)y + ---+{x:it))\ 

In section 2.c), we explain in details why we choose this drift. In this paper, we 
only prove a fluctuation theorem for the Gaussian case of this model: 

Theorem 1. Let > 0. Assume that 

Vx G R (p{x)=—-—T 
4(7^ 

and that, for any n > 1, the random variables X;^(0),..., Xn(0) are independent 
with common distribution p = We denote {lL{t), t > 0) the unique 

strong solution of the stochastic differential equation 

z^(t^ 

dz{t) =--^^dt + dB{t), z(0) = 0, (5a) 

where {B{t), t > 0) is a standard Brownian motion. Then, for any T > 0, 


/ SniVnt) 


, 0 < t < T 


n—¥-\-oo 


{uif), 0 < t < r), 


in the sense of the convergence in distribution on C([0,r],]R). 

This theorem suggests that, at least in the Gaussian case, our dynamical model 
exhibits self-organized criticality. Indeed it does not depend on any external 
parameter and the fluctuations of 5„( •) are critical: the processes evolve in a 
time scale of order and in a space scale of and the limiting process 
is the solution of the "critical" stochastic differential equation (Sa). This is 
the same behaviour as in the critical case of the mean-field model studied by 
Donald A. Dawson in [7], see section 3.a) for more details. 

For any n > 1, we introduce S* = + ■ ■ ■ + fff where • ,Cn) has the 

density proportional to 


(Xi, . . .,Xn) 


exp 


/ 1 (Xi H-h Xn)'^ 

y 2 x\ 4- * • * 4- 4- 1 


x|_+— 

2a2 ) 


with respect to the Lebesgue measure on R". In this paper, we also prove the 
following commutative diagram of convergences in distribution on R: 


SniVnt) 


(-4i) 
i —y ~hcx3 


j^3/4 


3 


(A) 


4- 

8 


3 


4- 

8 


(-42) 


U{t) 


t —y “hoo 

(-43) 



a 



-1 

exp 



ds 
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In section 2, we present some results on the general case of the model and 
we prove the convergences in distribution associated to the arrows (-4i), (-42) 
and (-ds) in the previous diagram. Next, in section 3, we give the strategy for 
proving a fluctuation result for our model and we explain that the Gaussian 
case is special because it can be analyzed through a two-dimensional problem. 
Finally we prove theorem 1 in section 4, i.e., the convergence in distribution 
associated to the arrow (-44). 

Acknowledgments. The author would like to thank the anonymous referee for 
his careful review and useful comments which helped to improve the presentation 
of the paper. 


2 Results on the general case of the model 

In this section, we first give general results on Langevin diffusions. Next we 
apply these results to prove existence and uniqueness of the solution of (Sa) 
and (S)?). We also prove the convergences in distribution associated to the 
arrows (-4i) and (-da). Finally we give a fluctuation theorem for an alternative 
version of the Curie-Weiss model of SOC. 


a) Langevin diffnsions 

Let / be a probability density function on M", n > I. The Langevin diffusion 
associated to / is a stochastic process which is constructed so that, in continuous 
time, under suitable regularity conditions, it converges to f{x) dx, its unique 
invariant distribution. 


Theorem 2. Let f be a positive probability density function on K.”, n > I, such 
that In/ is . We suppose that there exists K > 0 such that 

\/x G R” (V Infix), x) < K{1 + llxf). 


If (Bit), t > 0) is a standard n-dimensional Brownian motion and if f is a 
random variable in R" satisfying E(||/|p) < -boo, then there exists a unique 
strong solution to the stochastic differential equation 

dYit) = \v\nfiYit)) + dBif), iSf) 

with initial condition T(0) = /. Moreover {Y{t), t > 0) is a Markov diffusion 
process on R" admitting f{x) dx as unique invariant distribution and 


\/x G R” 


lim sup 


F{Y{t) G A|r(0) =x) 


f{z) dz 


= 0 . 


Proof. Theorems 3.7 and 3.II of chapter 5 of [9] imply that there exists a 
unique strong solution to (5/) with initial condition /, that its sample path is 
continuous and that it is a solution of the martingale problem for where 


\/g G (^^(R") 


^f9 


^/I5(ln/)\ dg 
(2 dx^ J 'bhi' 

t — 1 ^ I —I ^ ^ 
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Next, theorems 4.1 and 4.2 of chapter 4 of [9] imply that it is a Markov process 
and that its generator is {Af, D{Af )) with C D{Af). Finally theo¬ 

rem 2.1 of [16] gives us the uniqueness of the invariant distribution and the 
total variation norm convergence. □ 

Notice that this theorem is true if we remove the hypothesis that ^ has a finite 
second order moment, but the solution to (Sj) would be weak (see theorem 3.10 
of chapter 5 of [9]). 


b) Solution of (iS(j) 

Theorem 2 implies that (iSo-) admits a unique strong solution (U{t), t > 0) which 
is a Markov process whose unique invariant distribution is 


V2 


r 


a 




ds. 


Moreover 


lim sup 
t^+oo agBm 


F{uit)GA)-y^r(l 


exp 




= 0 . 


This is the convergence in distribution associated to the arrow (^ 3 ) in the 
diagram on page 3. 


c) Solution of 

In this subsection, we prove that (S[?) has a unique strong solution and that 
the convergence in distribution associated to (-4i) is true. 

Let us define „, the probability measure with density 


fn,p ■ y G 


1 / 1 (yi 

7 ^ exp 


+ Vn) 


2yi + --- +yi + l 




( 1 ) 


i=l 


with respect to the Lebesgue measure on M”, where Z* is a normalization con¬ 
stant. Let us prove that (S^) admits a unique solution. For any y G K", we 
denote 

Sn[y] =yi + - \-yu, T^y] =yl + -h 

and we notice that, for any j G { 1 ,..., n}. 


d 

dyj 


(1 iSn[y]f 

1 ^ 2 T„[y] + l 




2=1 


Sn[y] _ f Sn[y] 
Tn[y] + 1 ^^\T4y] + l 


2 

+ 2ip\yj). 


Therefore the system (S^) can be rewritten 

dXnit) = + dB{t), 

where B = (i?i,..., Bn). As a consequence, the solution of (S^) (if it exists) is 
the Langevin diffusion associated to f* p. 
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Let us introduce the operator L„ on (7^(18.") such that, for any / G C^(]R") and 
y G M”, 


Lnfiy) 


^ \ " ( 1 Sn [y] 

2 ^ j^\‘^T4y] + l 


^ ( Sn[y] 

2 Vr„[j/] + l 



dfjy) 

dVj 


Theorem 3. For any n> 1, there exists a unique strong solution 

(7c„(t), t>o) = ((xi(t),...,x;:(t)), t>o) 


to the system (S^) with initial condition X„( 0 ) having a finite second moment. 
Moreover it is a Markov diffusion process on K." with infinitesimal generator 
{L„, D{Ln)), where (^“(IR") C D{Ln), and whose unique invariant distribution 
is ffif p. Finally 

Vx G M" lim sup I V(Xn(t) G A I X„(0) = x) — ffif JA) I = 0. 

t^+oo I 


If we take p{x) = —x^jifia^) for any x G K., then theorem 3 proves the conver¬ 
gence in distribution associated to the arrow (Tli) in the diagram on page 3. 

Proof. Let n > 1. By hypothesis, there exists C > 0 such that 
Vx G M xip'{x) < (7(1 -I- x^). 


Moreover (p is (7^ on R thus the function In/* ^ is (7^ on R". For any x G R", 
we have 


r„[x](5„[x])^ ^ 

TnW+l {Tn[x]+iy ^ 


< 


(5'„[x])2 


+ 2C(n+||xf). 


(T„[x] + 1)2 

Next the convexity of 1 1 —>■ t"^ on R implies that 


VyGR” 


jSnly])^ nTn[y] 

(r„[y] + l)2 - (T„M + 1)2 - ■ 


since T„[-] < (T„[-] + 1)2. Therefore ff p satisfies the hypothesis of theorem 2 
and theorem 3 follows. □ 


Remark: we have chosen to built our dynamical model so that fif p is its unique 
invariant distribution. It is an alternative version of the Curie-Weiss model we 
designed in [4], given by the distribution 


dpnAxi,. ..,x, 


1 (l{xi 

0 = — exp ( - 


+ Xnf 


2 xf -I-h x2 


-f 2 >p{xi) dxi ■ ■ ■ dxn 


i=l 
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where Zn is a normalization constant. If we want to built the Langevin diffusion 
associated to the density of we obtain the system of stochastic differential 
equations 


dt + dBj{t) + i - Xl{t) ^ dt, 

j G {l,...,n}. 

In this case, the interaction function is not Lipschitz and we have to check first 
that Tn{t) ^ 0 for any t > 0: this would create technical difficulties to prove 
existence and uniqueness of a solution. In the next section, we give some results 
on the alternative version of the Curie-Weiss model of SOC (the model defined 
by the probability measure ^ - see formula (1)). 


d) The alternative Curie-Weiss model of SOC 

Let p be a probability measure on R. We consider an infinite triangular array 
of real-valued random variables (^^)i<fe<n such that for all n > 1, 
has the distribution 

dfi„Jxi,...,Xn) = ^exp [2 xl + --- + xl + l ) 

where Z* in the normalization constant. We define S* = -I- • • • -I- 

We obtain the same fluctuation theorem as in [11]. We only present the case 
where p has a density: 

Theorem 4. Let p be a probability measure having an even density with respect 
to the Lebesgue measure on R. and such that 

duo >0 / e^°^ dp{z) < -boo. 

dR 


If <j‘^ denotes the variance of p and pi its fourth moment then, under , 


^ A 

n—>-oo 


UaS j 


1/4 


ri- 


exp 


12cr8 ) 


ds. 


The proof of this theorem is given in section 18.b) of [12]. It is an adaptation 
of the proof of theorem 1 of [11], which consists in replacing the function F by 
the function {x,y) i—>■ x^ jifly -b 2/n). 

If we take '-pix) = —x^/(4ct^) for any x G R, then theorem 4 implies the conver¬ 
gence in distribution associated to the arrow (-42) in the diagram on page 3. 

3 Strategy of proof 

In this section, we first explain that the main ingredient for proving a fluctuation 
theorem for our dynamical model (in the case of a general function) will be 
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the study of its associated empirical process. Next we will focus only on the 
Gaussian case, i.e., when Lp : x i— > —x^/(4(7^) for some cr^ > 0. Indeed we will 
see that the Gaussian case can be handled by studying the convergence of the 
process 

/ / SniVnt) 

\ \ ’ 

We compute the generator of this process in subsection b). Finally we give the 
sketch of proof of theorem 1 in subsection c). 

a) The empirical process 

Let (p be such that has a unique strong solution ((X;^(t),..., t > O) . 

As in the equilibrium case (i.e., the alternative Curie-Weiss model defined in 
formula (1) or (2)), we would like to study the process (S„,Tn). However it is 
not Markov a priori, contrary to the empirical measure process Mn- It is the 
process taking its values on AIi(IR) and defined by 

-i n 1 ^ 

Vt>0 VAgBm = = 

‘ fe=i * fc=i 

where ((X^(t),... ,X^(t)), t > O) is the unique solution of 

Lemma 5. If the distribution of Xn{0) is invariant under permutation of coor¬ 
dinates, then {Mn{t,-),t > 0) is a Markov diffusion process on AIi(M). 

This lemma has a similar proof than lemma 2.3.1 of the article [7] - a paper by 
Donald A. Dawson about a mean-field model of cooperative behaviour. Daw¬ 
son’s model is defined through a Markov process which is solution of a system of 
stochastic differential equations. This process depends on two parameters and 
Dawson proves the existence of a critical curve in the space of the parameters. 
The critical fluctuations of the empirical measure process M„(-) evolve in a 
time scale of order and in a space scale of order We believe that our 

dynamical model has the same asymptotic behavior for the following reasons: 

* The invariant distribution of Dawson’s process is a particular case of the law 
of the generalized Ising Gurie-Weiss model, defined in [8]. 

* The alternative Gurie-Weiss model of SOG, defined in formula (1) or (2), has 
the same asymptotic behavior as the critical generalized Ising Gurie-Weiss model 
(see theorem 4). 

* The invariant distribution of our dynamical model is the law of the alternative 
Gurie-Weiss model (see theorem 3). 

Let n > 1. As in Dawson’s paper, we define the process by 

Vt > 0 VA G Br Un{t, A) = (^n(.\/nt, ^) ~ J dp{x) 

It takes its values on AI^(IR), the space of signed measures on R. 

The convergence of a sequence of Markov processes can be proved through 
the convergence of the sequence of their generators. Let us denote by the 


n 


1/4 


TniV^t) 


— a 


t > 0 



infinitesimal generator of Un- Let / and <I) belong to C^(R). We assume that $ 
is p-integrable. We have 


Vt>0 G^f (^J^<^{z)Un{t,dz)^ = 


where 


I X G 




If $ : z I—>• z then, for any i £ {1,..., n} and a; G R 


and 


dx^ 




If $ : z I—>■ z^ then, for any i £ {1,..., n} and x G 


dFf^^ 2xi d‘^Ff,<s> f N 4a:? 2 

and -^(x) = 


In both cases, if we suppose that (p : z i—>■ —z^/(4fT^), then we notice that, 
for any x £ M", the term LnFf^i^{x) only depends on n, 5'„[a:] and Tn[x\. This 
suggests that, in the Gaussian case, in order to prove the convergence of the 
process {^Sn{\/nt)/n^/^, t > O), we can turn the study of Un (which is a prob¬ 
lem in infinite dimensions) into a problem in only two dimensions. Indeed, we 
introduce the processes 5'„ and T„ defined by 


Vf>0 Sn{t) = = f zUnit,dz) 

Jr 

and 

Vt > 0 Tnit) = = f z^Un(t,dz). 

\ n J Jr 

In the rest of the paper, we suppose that ip{x) = for any a; G R. 


b) Generator of {Sn,Tn) in the Ganssian case 

Let n > I and f £ (R^). Let us define 4'/ on R" by 

Vx€R" = 

Proposition 6. For any n> 1 and f £ C^(R^), we have 

Vt>0 VJiLn'ffiXlit),...,X:i{t))=Gnf{Sn{t),fn{t)), 

where, for any (x,y) £ R^, 




dy 


2ct‘* dx 


2^ - 
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with 


Ri.{x,y) = y^{x,y)R^^\x,y) + ^{x,y) R^^\x,y) 


2x 52 / , ^ , 2y d^f ^ , 


9x9y ’ 9y^ ’ 

where (-Ri^^)ra>i and {R\?'^)n>i are sequences of functions from M? to K. verifying 

V/c > 0 lim sup max ( \Rl^\x,y) \, y)| ) = 0. 

”^+°° (x.y)GR" ^ ^ 

ll(a:,y)ll<fe 

Proof. Let us define '!'/ on K." by 

VccG R” 




Sn\x] Tjx] 


? 3/4 


- 


Let X £ R”. For any i £ {!,..., n}, we have 
d’I'f(x) ^ 


dyj 


1 (...) a/ 


a^d>/(a;) 1 a^/ 

2 ^3/2 1 ^ 3,2 


ay, 

where we write 


(...) 


dx ^ dy 

4x, ay 4x| ay 

^3/2 |9,;|9y ^3/2 ,5^2 

5„[a;l T„[a;l 


(•••), 




instead of 


^ 3/4 ’ 7j3/4 

in order to simplify the notations. We have then 


- ni/V2 


1 


Ln^f{x) = 


d^-^f{x) 


i=i 

1 av 


aa:2 


»S'n[a:] 
iPnia;] + 1 


— X, 


/ Sn[x] y a:A a^/(a;) 

Vr„[a;] + ly 0-2 I aa;, 


25„[a:] ay 2T„M av 


2-^71 (9x^ dxdy 

1 / n^/^S'„[a;] 


, 3/2 5y2 

-S'nN ^ 


2 \1 + T„[a:] n3/^(l + T„[a;])^ ) dx 

1/4 , 52 [a;] T„[a;] ^ a/ 


'^3/4(1 ^3/4^2y 


(•••)• 


We obtain that 


iTi i' ^ f ( ^"[2^] 1/4 2^ 

V^L„vI-/(x) = G„/ ^ -n/ a j, 

where Gnf is defined on R^ by 


y{x,y) £ 


Gnf{x,y) = 


2x a^/ 


dxdy 


2y 






2 aa;2 


{x,y) + ( -^^(1 - hn{y)) - ^hn(y) 


2(72 


2cr4 


aa; 


ay2 

(a:,y) 


v^y 


j,3/4jj4 


hn{y)'‘ 


ay 


(a;,y), 
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with 


hn -y & ]- +oo[i—;> + 


1 


na^ 


y 


We introduce the functions ^ and Sn' such that 


-( 2 ) 


Vy > /i„(y) = 1-+ -^el^\y) 

n^/^CT^ yjna^ yjn 


and Sn'^ (y) = hn(yY — 1- We obtain the formula of G„/ given in the proposition 
with 

Rn^-{x,y)' — > l^Sn\y) - ^£n\y) and : {x, y) ^ ^ ' 

It is easy to see that (i?n ^)n>i and {Rn^)n>i are sequences of functions which 
converge to 0 uniformly over any compact set in □ 


c) Sketch of proof of theorem 1 

Let us denote by the infinitesimal generator of the Markov process which is 
solution of (5 ct). It is defined by 

V/ G G2(R) Vx G M G./(x) = i/"(x) - 

Let n > I and f G (R.). By abuse of notation, we also write / for the function 
(x,y) G I —> f{x)- The essential ingredient for the proof of theorem I is the 
introduction of a suitable martingale problem. By Ito’s formula (see [17]), we 
prove that 

f(Sn{t)) = f{SniO)) + f G„/(§„(S)) ds + 

Jo 

where A^n,/ is ^ local martingale. By proposition 6, we have 

GnfiSn) = f'{Sn) + G./(§„)+/'(5„) i?« (5„, T^) 

'-V-" 

where (i?i^^)n>i is a sequence of functions which converges to 0 uniformly over 
any compact set in R^. 

Step 1: We notice that the term yl/(S'„,r„) does not converge a priori. To 
solve this problem, we introduce a perturbation: we transform the function / 
into a function Fnj which converges to / as n goes to oo, and which satisfies 

GjiFjif {^S — G (j f {^S “t“ a remainder. 

Notice that the perturbation theory and methodology was first introduced in [15]. 

Step 2 : For any A: > 1, we define the stopping time as the first exit time of a 
path of (S'n, Tn) from the domain [—k, fc]^, and we prove that j = Mnj{ - /\t^ /\T) 
is a martingale which is bounded over L^, for any T > 0 and fc > 1. 
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Step 3: We prove that P(r^ < T), the probability that a path of (^Sn, Tn) exits 
[—k,kY before the time T, goes to 0 when n and k goes to +oo. We also use 
the concept of collapsii^ processes (see appendix) in order to prove that the 
sequence of processes (T„(t), t > converges to 0 in the following sense: 


Vry > 0 


lim P( sup |T„(t)| > 77 ) = 0. 

n-s-+oo \o<i<T / 


Step 4: We prove that the sequence (5'„(t), t > 0)„>i is tight in the Skorokhod 
space T], R). 

Step 5: We deduce from the previous steps that there exists a subsequence 
which converges in distribution to some process U on 2?([0,T],R). 
We prove then that, for any fc > 1 and t G [0, T], 






n—¥-\-oc 


ptAT 

Mf{t)= f{U{t A T)) - f{um - / G„f{U{s)) ds 

Jo 


and that At/ is a martingale. As a consequence U is uniquely determined as 
the unique solution of the martingale prcAlem associated to Ga. We conclude 
that U is the solution of (5^) and that converges in distribution to U 

on 2?([0,T],R), and thus on 

These steps are developed in detail in the next section. 


4 Proof of theorem 1 

Step 1: Perturbation 

Let f G G^ (R). We want to find functions H / and K / defined on such that 
Kj ■■ {x,y) I—^ fix) + -^Hf{x,y) + -^Kfix,y), 

Ti ' y Tl 

satisfies _ _ 

GnFnJ = Gaf + Rnji 

where Rn.f is a remainder term. Let us find necessary conditions. We suppose 
that we have built Hf and Kj and we assume that they are We have then, 
for any [x, y) G R^, 


GnFnj{x,y) = n 


= r,l/4 


y dHf^ 

0-2 dy 


ix,y) - -^fix) - 


y dKf 
(T^ dy 


ix,y) 


^ ~ + i/"(x) + a remainder. 


2ct‘* dx 

The function iJ/ should verify 


V(x,y) G 




dy 


2cr4“ 


We choose 


Hf : {x,y) I—^ 
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Therefore the function Kf should satisfy, for all {x,y) G 


GnFnj(x,y) = 

+ [xy^ — x^CT^) f'{x) + ^f"{x) + the remainder 
^ ^ ^(3/'(a;) +a:/"(x)) 

a;3 

~ TT^f'i^) + + the remainder . 

2(7^ 2 

So that the variable y disappears in the leading term of GnFnj(x,y), the func¬ 
tion Kf should verify 


V(x,y) G + |^(3/'(x) +xf"{x)) = 0. 

We choose 

Kf : {x,y) I—^ |^(3/'(x) -f x/"(x)). 

It is easy to see that these choices iorjl / and Kf are sufficient for the variable y 
to disappear in the leading term of GnFnj{x,y). The remainder term is then 


Rnj — Rn 


ri/4 


K 


Hf 


\/n 


K 


Kf 


We notice that, so that the above computations are possible, it is necessary that 
/ is C^. Indeed, the hrst four derivatives of / appear in the remainder term. 
We also remark that, if / G C"*(]R), then the functions Hf, Kf and their hrst 
and second derivatives are bounded over any compact set in Finally let us 
recall that (i?n^)n>i and {Rn^)n>i are sequences of functions which converge 
to 0 when n goes to -boo, uniformly over any compact set. As a consequence we 
have the following proposition: 

Proposition 7. Let n> I and f G C''^(]R). We define Hf and Kf o?z by 


V(x,2/) G 


Hf{x,y) = 


Kf{x,y) = |^(3/'(x) +xf"{x)). 


Then the function 


Fnj ■■ {x, y) 


fix) 


^Hf{x,y) + -^Kf{x,y), 


verifies GnFnj = G^f + Rnj, with Rnj o, remainder term satisfying 


\/k>0 


lim 

n—>-+oo 


sup 

>,y)eR^ 

(x,v)\\<k 


Rnjix,y) 


= 0 . 
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Step 2: Introduction of a martingale problem 

We give ourselves n > 1 and / G For any t > 0, we have 

/ = /(»«(')) = (sM.f,.{t)). 

We define the process i > 0) by 


vt > 0 Mnjit) = Fnj{Snit),T„{t)) - F„j(0), T„(0)) 

-[ GnFnj{Snis),fn{s))ds. 
Jo 

By applying Ito’s formula to the function 


^ n,f • ) G 


F 


Sn [^] Fn [x] 2 


I n3/4 ’ n3/4 


we obtain 


^ {•'t 

Vt > 0 Mnjit)=n^^^^ (X„(^/ns)) dBj{s). 

„-_i ^0 

It is a local martingale and 


i=i 


n pt 


Vt > 0 {MnJ,Mn,f)t = 


f=l ' 




n,/ 


dxj 


(X„(i/ns)) ds. 


For any fc > 0, we introduce the stopping time defined by 
= inf I |5„(t)| > fc or |T„(t)|>fc|. 

Let r > 0. We denote j{t) = At^ AT) for any t > 0. 

Lemma 8. For all k > 1, n > 1 and f G C3(R), the process *3 a 

martingale which is bounded over . Moreover 

Vt > 0 sup E f{t)^) < + 00 . 

n>l 


Proof. For any t > 0, we have 


{M 


k 

nj 





2 

(XniVns)) ds. 


Moreover, for alH G {1,..., n} and x G R.", 


, _ 1 r/ ( 

dxi 71 - 3/4 71 , 3/4 

1 / 1 dHf 1 dKf\ 

17.3/4 \^j^l/4 ^1/2 j 


>5'n[x] r„[x] 
7i3/4 ’ 7i3/4 


- ni/V" 


2x, / 1 diL/ 
7i 3/4 \^n3/4 dy 


1 dKf \ ( 5'n[x] Tn[x\ ^ 1/4 2 

ni/2 dy ) \ n3/4 ’ n3/4 


( 3 ) 
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By squaring these terms and by summing over all t G {1, ... ,n}, we observe 
that there exists a constant C* > 0 such that, for all x G M" verifying 


■5'nN 

^3/4 


< k 


and 


Tn[x] 

74 , 3/4 


- n^/V^ 


< k, 


we have 



As a consequence, for any t > 0, 

snp E {{Mlf,Mtj)t)<C^T. 

n>l 

Therefore, for any n > 1, the process M^ ^ is a martingale bounded over 
(see theorem 4.8 of [14]) and 

Vt > 0 E {Mtjitf) = E {{Mlf, Mtj)t) < C) T. 

This ends the proof of the lemma. □ 


Step 3: Study of the asymptotic behavior (r,^)„>i 
Lemma 9. For any £ > 0, there exist rig > 1 and fcg > 1 such that 

sup P {t^'^ <T) <s. 

n>ne 


Moreover the process {Tn{t), t > collapses to zero, i.e., 

Vry > 0 lim P( sup |r„(t)| > ry ) = 0. 

n^+oo \o<t<T ) 

Proof. Let k,e > 0 and n > 1. We have 

p(T„"<r)<p( sup |f„(t)|>M+p( sup |5„(t)|>M. 

Y0<i<TAT^ ^ / \0<t<TATyj ^ / 

We denote P(A*) + P(B*) the sum in the right side of this inequality. 

Let us deal with the bound of P(Aj)). To this end we would like to apply 
proposition A.2 in appendix to the positive semimartingale {^n{t), t > 0)„>i 
defined by ^ 

Vn>l Vt>0 ^nit)=Tn{tf. 

By applying Ito’s formula, we get 

d^n (t) = G„/o (§„ it ), f„ (t)) dt + Tn (t) dB, (t ), 
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with /o j {x,y) I— >^y^- With the notations of proposhion A.2, we have 
Cn(0 = Gn fo(Sn(t),Tn{t)) and Zn,i{t) = 4n“(Vnt) T„(t) for all t > 0, 
n > 1 and z G rz}. We have 

n 1 ^ 

Vn>l VtG[0,r^] = 

Hence condition (‘^ 4 ) of proposition A.2 is verified with = 16fc^(CT^ + k). 
Next, by proposition 6 , for any n > 1 and t G [0, 

Ut) = +4^2+2 r„(t)i?i 2 ) (A„(t), fnit)) + - 47 jTn{t) 

<+40-2+ 2 fc sup R^^\x,y) 

o- \\{x,v)\\<k n / 

Condition {'Z’ 3 ) is then verified with k„ = ^/ri for any n > 1, C 2 = ‘2/a^, 

(74 = 4ct^ + 2k sup sup R^\x, y) + 4fc < +00 
">i ll(a:,y)ll<fe 

and C 3 , (/3n)n>i may be chosen arbitrarily. We choose (/3n)n>i such that fin j kin 
goes to 0 when n goes to + 00 . 

Let us examine condition ('^ 2 ): we denote Yf = (AT* (0))^ — cr^ for any 
z G {1,..., zz}. Since Ar;^(0),..., Ar”(0) are independent random variables with 
common distribution A/’(0, cr^), we get that Yff are independent identi¬ 

cally distributed random variables which are centered and have finite moments 
of all orders. Theorem 2 of [3] implies that, for any v > 2, there exists > 0 
such that 

Vrz>l E(|y„i+ --- + y”|") < iG,rz"/2. 

Hence, for all d > 1 and rz > 1, 

e[(«..(0))''] =E + + = 

Condition ('^ 2 ) is then satisfied for any d > 1, with Ci = K 2 d and a„ < ^/ri for 
all rz > 1. So that condition (‘^ 1 ) is verified, we choose d > 2 and = rz^/^ for 
all rz > 1. We have 

V antif^ = V rz'^/^ = 

As a consequence, proposition A.2 implies that there exist M > 0 and rzi > 1 
such that 

suppf sup \fnit)f>MnG(^‘^^-G^\<i, (4) 

n>ni \0<t<TAT^ I ^ 

We increase the value of rzi so that 
sup P sup 

n>ni \0<i<TAT^ 
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Let us deal now with the term In the rest of this proof, we assume that 

/ is the function {x,y) i—>■ x^. We have 

W ^ 1 O 17 ta l4.\ 'L' I Tn{t) Sni^Y'Tnit)"^ 

Vn>l Sn{t) =Fnj[Sn{t),Tn{t))+—-^^ - -^=^3 -, 

thus 

Vn > 1 F^j{Sr.{t),f^{t)) = S^itY (l - 

We obtain that, for n large enough, 


F{B^)=f( sup 

\ 0<t<TAT^ 4 


< 


<P sup Fnj{Sn{t),Tnit)) > — ] 

\ 0<t<TAT* o J 

Fn,f{Sni0),Tn{0)) > l +P| SUp AinjY) > 


+ Pj sup GnFnj{Sn{t),Tn{tY > 

\ 0<t<TAT^ 


For any n > 1, the random variables X^{0),..., X^{0) are independent with 
common distribution ^^(0, cr^) thus, by the Central Limit Theorem, we get 
(S'n(0))„>i and (T„(0))„>i converge in distribution to 0. This implies that, 
for n large enough, 


F„,y(^„(0),T„(0)) >- 



Next proposition 7 gives us 


GnFnj{Sn{t),Tr,{t)) =1 

< 1 


\Rnj{Sn{t),f4t))\. 


and 


lim sup \Rn f{u,v)\ = 0. 

\\{u,v)\\<k 

If we choose k > '/2YT and n large enough, then 


sup GnFnj{Snit),Tn{t)) > 

0<t<TAT* 


< p 


1+ sup \Rnj{u,v)\> — 
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Finally, by lemma 8 , j 
implies 

P I sup 

Y 0<i<TAT^ 


is a martingale thus Doob’s maximal inequality 





- (fcV24)2 


Lemma 8 also implies that (E(A1 ^j(T)^))^ is a bounded sequence. Hence, 
for k large enough, 


P [ sup 

Y 0<t<TAr^ 






As a consequence, there exist n 2 > 1 and kg > 1 such that P(i3*®) < e /2 for all 
n > 712- We denote Ug = ni V n 2 . We have proved that 

Vn > ng P (rt <T)< P(A^= ) + P(H^') < e. 

Let us prove the second assertion of the lemma: for any 77 > 0, we have 


P( sup |rn(f)| > ?7 ) < P ( sup |T„(t)|% ryM +P(t^' < r) . 

\0<t<T J \o<t<T/\Tn^ J 

By formula (4), for n large enough, 

pf sup |f„(f)|>, 7 ) <^+P(r^<T)<^. 

By letting e goes to 0, we obtain that (r„(f), t > collapses to zero. This 

ends the proof of the lemma. ~ □ 

Step 4: Tightness of {Sn{t), t > 0)n>i in 2?([0,T],IR) 

Since {Xn{t), 0 <t < T), n > 1, and the limiting process 0 < t < T) be¬ 

long to (^([0, T], R), it is enough to prove that (5'„(t), t > 0)„>i is relatively com¬ 
pact for the weak convergence in ^([O, T], K), which is a Polish space (see theo¬ 
rem 12.2 of [2]). Prohorov theorem (theorem 5.1 of [2]) implies that it is enough 
to prove that (5'„(t), t > 0)„>i is a tight sequence. As in [6] and [5], we use the 
following tightness criterion: 

Proposition 10. A sequence 0 < t < T)„>i on 21([0,T],M) is tight if 

(a) for any s > 0, there exists M > 0 such that 

sup P [ sup |Cn(i)| > m \ < e, 

n>l \0<t<T ) 

(b) for any e > 0 and rj > 0, there exists <5 > 0 such that 

sup sup P ( |Cn(T' 2 ) - Cn(n)| > ?7) < 

n>l ri,T 2 G%i 

0<Ti <T2 <(ti+(5)AT 

where, for any n > 1, Tn is the set of all the stopping times adapted to the 
filtration generated by the process 
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Lemma 11. The sequence (Suit), 0 <t < is relatively compact for the 

weak convergence on 21([0,T],R). 

Proof. It is enough to prove that {Sn{t), 0 < t < r)„>i verifies conditions (a) 
and {b) of proposition 10. In the proof of lemma 9, we proved that, for any 
a > 0, there exists ka > 0 and Uq > 1 such that 

sup P <T) < a 

n>na 

and, for all n> Ua, 

P { sup |5„(t)| > 

\0<t<TAT^“ ^ ^ 

We give ourselves e > 0 and we denote a = 2e/?>. We obtain that, for all n > Uq,, 


P ( sup \Snit)\ > % ) 

\0<t<T ^ J 

<f( sup |5„(t)| > +P(t„^“ <T) < ^=£. 

\0<i<TAT^® ^ ^ 

Hence condition (a) is verified. 

We prove now condition (6): we give ourselves n > 1 and e,ri,S > 0. Let 
Ti and T 2 be two stopping times adapted to the filtration generated by the 
process 5'„ and such that 0 < ri < r 2 < (n + d) A T. Setting a = 2e/3, we have 

P(| 5 '„(t 2 ) - 5'„(ri)| > rj) 

< P (|5„(t 2 A r^) - §„(ti a r^) I > 77 ) + P (r^ < T) 

< - E ( |5'„(t2 a r*“) - Sn{Ti A r*“)|) + a, 

where we used Markov’s inequality. In the rest of this proof, we assume that / 
is the function {x, y) 1 —> x. We have 


|5'„(t2 a - 5'„(ti a 

^ A r„"“),r„(r 2 A A A r„^“))| 

{Sn (t2 a ), f „(t2 a )) - Ky (§„ (n A ), f „ (n A r^)) | 

\GnFnj{Sniu),f^iu))\ du + 




pT2 Ar^ 
/Ti ATn“ 


We have 
E 


(Af^j(r2) 


= E E 
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where Qp = a(A4^y(s), 0 < s < t) for all t > 0. By lemma 8, is a mar¬ 

tingale bounded over thus it is uniformly integrable. Martingale Stopping 
Theorem (theorem 3.16 of [14]) implies that 


Hence 


Mty(ri)=E 





E 


(Miyiry-M 


nj 


= E 


Mtyir2f 


+ Mty{nr-2Mty{ri)E 


Miy{T2) 


= E 


M'pyi.r.f 


-Miyy.f 


and thus 

e/^e 




-m: 


= E E 


(Mty,Mtyy. - (Mty,Mty}n 


= E Vn 


r2 Ar 


7 = 1 ^"^lATn 


dF, 




dxi 


(^XniVnu)) duj <Ci°‘S, 


where Cj“ is the constant introduced in the proof of lemma 8 for fc = fco,. We get 


E 


(\M^yyT^)-Mly{Ty\)<^^5. 


Next, since / : {x,y) i —> x, proposition 7 yields 

Snity 


and 

Therefore 

rT'2r\Tp° 

I ^ ka, 

TiAr^i 

Finally 


GnFn,f{Sn{t),Tnit)) = - ^ + RnJ (S n{t), T n{t)) 


V/c > 0 lim sup \Rn,f{x,y) \ = 0. 

"^+°° Wix,y)\\<k 


\GnFr^j{Sniu),Tn{u))\ du < ( ^ + sup \R„j{x,y)\\S. 

\\G,v)\\<ko. J 


\Hf{Sr.{T2 A T„^“),T„(r2 A T^y) - Hf{syn A r„^“),T„(Ti A | < 


and 


Kf{SyT2 A r„"“),T„(r2 A t^)) - Kf{SyT, A r„"“),T„(Ti A T^)) | < 


dcr^ 
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Hence, for n large enough and S small enough, 

E ( A - S^(n A r„^“)| ) < ^. 

We obtain 

P ( \Sn{T2) - 5„(ti)| >r])<^=e. 

Condition (b) of proposition 10 is then satisfied and this ends the proof of the 
lemma. □ 

Step 5: Identification of the limiting process and convergence 

Let us identify the limiting process. By lemma 11, there exists a subsequence 
t> which converges in distribution to some process (W(t), t > 0) 

on P([0,T],R). By lemma 9, t > converges in distribution to 

the null process on 2?([0, T], K). 

For A: > 0, we introduce the stopping time 

= min (^T, inf | |T„(t)| > ^}^ • 

If t > T then P(t,^ < t) = 1 and, if t < T, then 


lim P(t,^ <t)< lim P ( sup |r„(t)| > fc ) =0, 

ri—>+oo n^+oo y 0<t<T / 

by lemma 9. As a consequence (r,^)n>i converges in distribution to T. 

We give ourselves / G C''^(IR). For any n > 1 and t € [0, T], 

Pn,f{Snit),Tn{t)) = /(5'„(t)) + ^ {Snit) ,T n{t)) , 

the functions Hf and Kf being continuous. Next, proposition 7 implies that, 
for any n > 1 and t G [0, T], 

= G^f{Snit)) +Rnj{Sn{t),fnit)), 
where Rnj is a continuous function on such that 

V/c > 0 lim sup |i?ri f{x,y) \ = 0. 

\\G,y)\\<k 


Let fc > 0. For any t > 0, we obtain 

A Mf{t) = fmAT))-fiUi0)) 

n—>-+oo 

For all n > 1 and t G [0, T], we have 


rt/\T 


G^fiUis))ds. 




tAT„ 


1=1 ■ 






dxj 


(XniVns)) ds, 
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and, using formula (3), we get 


j=i 


= f/'(sj + 

'4T„ 


1 dHf 1 dKf 




+ 


dx 9a; 

1 9i?/ 1 dKf 


4^„ 


n 


1/4 


4(7 

/'( 5 „) 


(9y dy \ 

1 ai// 1 dKf 
v}!^ dx 71^/2 9a; 

1 9i7/ 1 dKf 


{Sn,T, 

{Sn, fn) 

(§„,T„) ) . 


9y dy 

Assume that / has a compact support. Then we observe that there exists a 
constant Cf such that 


Tn{t) < k 


v^i: 

/=i 


/ dlpnj 
V dXn 


(X„(v^t)) < Cf 


As a consequence Mnj{ - A rf) is a martingale and 

Vt>0 sup E(7W„j(t At,J)2) < <+ 00 . 


n>l 

This implies that, for all t > 0, A is an uniformly integrable 

family. Therefore A4/ is a martingale. 

Theorem 1.7 of chapter 8 of [9] implies that the martingale problem associated 
to { (/, Gaf ) ■ f G (^“(IR)} admits a unique solution: it is the strong solution 
of the differential stochastic equation 

dz{t) = — A dt + dB{t), 2 ( 0 ) = 0, 

2(7'* 

where t > 0) is a standard Brownian motion. As a consequence the 

limiting process 0 < t < T) is uniquely determined. Therefore 

f^.o<t<T) = (s„(t), 0 < t < r)„^, 

n / n>l 


converges in distribution to (9^(t), 0 < t < T) on T>([0,T],R). Finally, since 
the sample paths of {K{t), 0 < t < T) are continuous, this convergence in 
distribution holds in C'([0,T],R). This ends the proof of theorem 1. 


Appendix 

A proposition on collapsing processes 

Definition A.l. A sequence of real-valued stochastic processes {£,n{t), t > 0)„>i 
collapses to zero if 

Ve > 0 VT > 0 lim P ( sup |^n(t)| > e ) =0. 

n^+00 \o<t<T ) 
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The concept of collapsing processes has been developed by Francis Comets and 
Theodor Eisele in [6]. 

Proposition A. 2. Let {^n{t), t > 0)„>i be a sequence of positive semimartin¬ 
gales on a probability space P). For any n > 1, we give ourselves an 

integer mn > 1 and independent standard Brownian motions "which 

generate a filtration (-Ft)j>Q- We assume that there exist {Ft)adapted pro¬ 
cesses {Cn{t), t > 0) and {Zn,i{t), t > 0)i<i<mn such that 


dfn{t) = Cn{t)dt + ^ Zn,i{t)dBi{t). 
i=l 


We suppose that there exist d > 1, positive constants Ci,... ,C^, increasing 
sequences (k„)„>i, (a„)ra>i, {Pn)n>i (in-d a sequence (r„)„>i of stopping times 
verifying 


Una 


-1 

n 


>+oo 


Oin ^ 0, 

n—¥-\-oo 


-1 


Pn 


>+oo 


Vn > 1 

Vn > 1 Vt G [0, Tn] 



Ut) < - 


< Clan‘s, 

nnC2^n{t) + PnC^ + C4, 


and 

rrin 

Vn > 1 Vt G [0, Tn] ^ Zn,i{t)'^ < Cs- 
Then, for any ^ > 0 and T > 0, there exist M > 0 and no > 1 such that 


C^i) 

C^3) 

(^^4) 


sup 

n>no 


sup ^n{t) > M I 

0<t<T/\Tn 


a„ ^ V OriK, 


-1 


< e. 


This is proposition 4.2 of [5]. It is a simple adaptation of the proposition in 
appendix of [6]. 
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